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Abstract 

A simple and consistent real time analysis of the long-wavelength chiral conden- 
sate fields in the background of hard thermal modes is presented in the framework of 
43 the linear sigma model. Effective evolution equations are derived for the inhomoge- 

neous condensate fields coupled to a heat bath. The multiple effects of the thermal 
background on the disoriented chiral condensate are studied using linear response 
theory. We determine the temperature dependence of the equilibrium condensate, 



and examine the modification of the sigma and pion dispersion relations as these 
mesons traverse a hot medium. We calculate the widths by determining the dissipa- 
tive coefficients at nonzero temperature at one- and two-loop order with resummed 
meson masses. Our results show that not only decay processes, but elastic scattering 
processes are significant at high temperatures, yielding to short relaxation times in 
the phase transition region. The relaxation times obtained are shorter than in pre- 
vious estimates, making the observation of DCC signals questionable. Throughout 
this work Goldstone's Theorem is fulfilled when chiral symmetry is spontaneously 
broken. 
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1 Introduction 



Non-equilibrium phenomena in many-body systems has received a great deal of attention 
during the last two decades. The theoretical interest keeps growing as new experiments 
are readily emerging in different research areas of physics. Such are the chiral phase 
transition and quark-gluon plasma in ultarrelativistic heavy ion collisions, ultrafast spec- 
troscopy in semiconductors, Bose-Einstein condensates, strongly correlated Fermi systems 
in condensed matter physics, and electroweak baryogenesis and inflation in Early Universe 
cosmology. In the following we present a quantum field theoretical description of the dy- 
namics of nuclear matter formed as a consequence of nuclear collisions at ultrarelativistic 
energies. 

The existence of a deconfined, chirally symmetric phase of QCD was predicted long 
ago [p]]. The interplay of results obtained from the three major research approaches, 
lattice simulations, theoretical models and experiments at RHIC-BNL, lead us to expect 
that 'finding' quark-gluon plasma is within reach. PHENIX results || from the Year-1 
run of RHIC suggest that the energy density achieved in Au+Au collisions already at 
yfs = 130 GeV is high enough to be favorable for the existence of free quarks and gluons. 
Early thermalisation of the partonic medium || has been indicated by the elliptic flow 
analyses at STAR [|J. However, most probably thermal equilibrium is not maintained as 
the matter rapidly expands and the temperature quickly drops through the QCD phase 
transition. Arguments that hint toward non-equilibrium, even explosive || dynamics in 
heavy ion collisions, such as large fluctuations of the average transverse momentum, and 
almost equal sideward and outward HBT radii, were reported in |J. The dynamical 
evolution of such an out-of-equilibrium system is not yet understood. When it comes 
to possible approaches, lattice simulations unfortunately do not prove to be the way to 
go: lattice can describe static situations in thermal equilibrium. Therefore, for learning 
about the dynamics of non-equilibrium systems, one has to reside to different theoretical 
approaches. Attempts are made using field theory |7], ||, |, |K], |Il|] , covariant kinetic theory 

and recently also non-equilibrium fluid dynamics f!3|| . 

Until now no unambiguous observable that could serve as evidence for the phase 
transition has been identified. Based on the knowledge that a transition from ordered 
to disordered phase is accompanied by the formation of condensates, disoriented chiral 
condensates (DCC) have been suggested |TJ| as signature of the chiral phase transition. 
This means that the matter formed in a heavy ion collision can relax into a vacuum 
state that is oriented differently than the normal ground state. Relaxing of DCCs to the 
correct ground state then happens through the emission of low momentum pions with an 
anomalous distribution in isospin space. Detecting fluctuations in the ratio of produced 
neutral pions compared to the charged ones can serve as a signal for the deconfinement and 
chiral symmetry restoring phase transition. Such DCC signals have not been observed at 



CERN-SPS energies [15|. The STAR detector at RHIC searches for dynamical fluctuations 



on an event-by-event basis and may be therefore better suited for DCC searches. Other 
signals were proposed in context of DCCs: dileptons | 16| , and recently the anomaly in 
the Omega and anti-Omega abundances observed at SPS ITT . 



The ability to detect DCCs depends on their lifetime. In the original work jTJJ DCC 
formation was proposed with the assumption of a perfect quench. This means that after 
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the critical temperature is reached the long-wavelength modes of the chiral condensate 
decouple from the thermal modes and evolve according to zero temperature equations 
of motion. A more realistic discussion, accounting for the presence of a thermalized 



background was first presented in |18| . Due to possible energy exchange between different 
degrees of freedom dissipation occurs, which in turn reduces the lifetime of the condensate. 
This problem has been addressed previously in |7|, ^ |^, |TD|, [II], [B5], |2D[. Most of the 



previous calculations though focus on the evolution of the order parameter only 0, [TT|1 . 
Based on the observation that true pions have a finite width ET| the damping of pion 



condensates was considered || pLOfl , but only at one- loop level. The usual argument for 
neglecting two-loop contributions is that these are higher order in the coupling. However, 
there is an important reason for why two-loop contributions are not negligible: whereas 
decay processes can happen only under some kinematic conditions, scattering can always 
happen. Dissipation due to scattering has been evaluated for the scalar condensate in 
"|Tl]| . A first attempt to look at the lifetime of DCCs based on scattering processes was 



presented in [T9[ . 

In this paper we present a quantum field theoretical analysis of homogeneous and 
inhomogeneous chiral condensate field configurations out of thermal equilibrium that are 
coupled to a thermal bath. The framework is the linear sigma model, which has proved to 
describe fairly well the hadronic phase of two-flavor QCD. We treat the long wavelength 
chiral fields classically and we account for the short wavelength fields perturbatively, but 
resumming certain diagrams. Such semi-classical description is acceptable since the occu- 
pation number of the low momentum modes is large. The effect of the heat bath on the 
condensate is contained within the deviations of the thermal field fluctuations from their 
equilibrium values. We identify these deviations as time-delayed responses of the hard 
thermal modes to the presence of the condensate, and evaluate them using linear response 
theory, which provides the evolution of observables in real physical time. The response 
functions renormalize the equations of motion, modifying the particle properties, and give 
rise to dissipation. Due to possible interaction and thus energy exchange between the soft 
and the hard modes decay channels open up and particles can scatter. These processes 
are responsible for the dissipation of the condensate. To the best of our knowledge, a 
consistent incorporation of relaxation processes at one- and two-loop level has not been 
done before in this context. There is another important effect, that of the change in the 



velocity. We discuss this in p2 



It is important to determine the order of the chiral transition, as this influences 
the dynamical evolution of the system. Experimentally, large-acceptance detectors are 
now able to measure average as well as event-by-event observables, which in principle can 
distinguish between scenarios with a first order, a second order, or merely a crossover type 
of phase transition. QCD with two massless quarks possesses chiral symmetry, described 
by the SU(2)r x SU{2)l ~ 0(4) group. This continuous symmetry is spontaneously 



broken at low temperatures, resulting in the emergence of Goldstone bosons ||23|| . Based 
on universality arguments p|J the transition between the chiral symmetric and symmetry 
broken phases is of second order. Nature, however, supplies a different situation, with 
non-zero quark masses. Small quark masses explicitly break chiral symmetry, which then 
alters the nature of the phase transition. The second order transition becomes a smooth 
cross-over, provided that the baryon density is zero. In case of a sufficiently large baryon 
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chemical potential the transition is of first order, suggesting the existence of a tricritical 
point in the T) plane of the QCD phase diagram [p5| . Year-1 RHIC results suggest 



that at -y/i = 130 GeV the baryon density is approaching the low density limit. Expecting 
an even smaller baryon density at the maximum collision energy justifies us to work at 
finite temperature and zero baryon chemical potential. 

The paper is structured as follows: In Section |2| we derive coarse-grained equations 
of motions for the long wavelength chiral condensate fields which are in in contact with 
a heat bath. In Section |3| the response of the thermal bath to the presence of the non- 
thermal condensate is evaluated. In Section |] we present the temperature-dependence of 
the equilibrium value of the condensate and that of the self-consistently calculated meson 
masses. It is important to emphasize that our analysis is performed with the fulfillment 
of Goldstone's Theorem for all temperatures below the critical one. Also, the tachyon 
problem appearing in the mean field treatments is eliminated. In Section |5| we analyze 
the dissipation of long-wavelength sigma and pion fields due to decay and scattering 
processes. In Section ^| the relaxation time of DCCs is determined. We summarize our 
results in the concluding Section [F[ 



2 Deriving Equations of Motion 

Within the linear sigma model framework the study of DCC formation and evolution is 
convenient, since the physical picture is rather transparent. The theory is formulated 
in terms of the chiral field $ = (cr, 7r). The scalar sigma field a describes the scalar 
quark condensate (qq), which serves as order parameter. The pseudoscalar pion field 
7T = (7T 1 , 7T 2 , vr^" 1 ) is directly related to the pseudoscalar condensate, (qf^q). The 
dynamics of the chiral condensate is completely determined by the evolution equations in 
space and time for the long-wavelength chiral fields. In the following we derive effective 
equations for these low-momentum fields in a hard-momentum background. 
The Lagrangian of the linear O(N) sigma model is 

= h^d^ - u{$) , (l) 

where the usual choice for the potential is 

U(<S>) = ±(&-vtf + Ha . (2) 

The explicit breaking of chiral symmetry is implemented through the H-term that tilts the 
potential in the sigma direction. H = f w m% and Vq = f%—ml/\, where f n = 93 MeV is the 
pion decay constant, A is a positive dimensionless coupling constant, and m n = 138 MeV 
is the zero temperature mass of the pion. The parameters of the Lagrangian are chosen 
such that for H = chiral symmetry is spontaneously broken in the vacuum. The 
potential then resembles the bottom of a wine bottle with the zero temperature minimum 
at $ = (/,r,0). The excitations in radial direction, the sigma mesons, have a mass 
ml = 2\f%, and excitations along the azimuthal direction, the pions, are Goldstone 
bosons, with = 0. 
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The evolution of the fields in thermal equilibrium has been extensively studied 
after the pioneering works by Linde |28[] and Dolan and Jackiw ||29|| . The usual procedure 
is to express the fields as 

a(t,x) = v + (jf(t,x) , 

n%x) = n)(t,x) , i = l,...,iV-l . (3) 

The thermal average of the chiral field, $ = (y,0), is the order parameter chosen to 
lie along the sigma direction, (cr) eq = v and (v?) e g = 0. The fluctuations about the order 
parameter average to zero, (07) eg = (^f)eq = 0. Throughout this work we use the following 
notation: 

• (O) is the non-equilibrium ensemble average of an operator O ; 

• {0) eq denotes the equilibrium, but interacting ensemble average; 

• (O)o is the free ensemble average. 

For non-equilibrium thermal conditions, the idea is that instead of writing the fields as 
in (H), we separate them into their low and high frequency modes. This can be done, for 
example, by introducing a momentum cutoff A c , as in f3T| . Then, by integrating out the 
high frequency modes we obtain an effective theory for the low frequency, long wavelength 
modes. The fields can be separated as follows: 

a(x) = a(x) + CTf{x) 

tx\x) = tP'O) +ir)(x) , i = l,...,N-l (4) 

• a and 7p are slowly varying condensate fields, representing low frequency modes 

— # 

with momentum | k |< A c . These soft modes are occupied by a large number of 
particles and may then be treated as classical fields; 

• <jf and 7Tj are high frequency, fast modes with | k |> A c . These hard modes, 
representing quantum and thermal fluctuations, constitute a heat bath. 

A choice of A c = would mean that only homogeneous condensates, with k = 
momentum are studied. Here we discuss condensates that can also be inhomogeneous. 

The problem to be solved, then, is to describe the evolution of long wavelength classical 
fields that are embedded into a thermal bath. In our approach, these soft modes follow 
classical equations of motion, whereas the effect of the hard thermal modes is taken into 
account in a perturbative manner. We are well aware of the problems of a perturbative 
treatment, since the sigma model is a strongly coupled effective theory: With the choice 
of for example m a ~ 600 MeV for the vacuum mass of the sigma meson the coupling 
constant is A ~ 20. However, we improve the model by resumming a certain class of 
diagrams in the perturbation series. We believe that even if the self-consistent solutions 
are approximate only, they still yield to qualitatively reliable results. 

We average the Euler-Lagrange field equations over time and length scales that are 
short compared to the scales characterizing the change in the slow fields, but long relative 
to the scales of the quantum fluctuations. This is known as coarse-graining. The average 
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of high frequency fluctuations is thus (<T/(x)) = and (iff(x)) = 0, while {v(x)) = a(x) 
and (tt 1 (x)) = tt 1 (x). It should be noted at this point that, compared to earlier works 



(for example |TT|), we allow for a nonzero ensemble average not only along the sigma 
direction, but also in the pion directions. In other words, we allow the formation of 
disoriented chiral condensates. Also, cross correlations between fluctuations of different 
fields are nonzero, {fifj) ^ 0, and moreover, we consider nonzero cubic fluctuations of 
the form ((pi(pj(pk) 7^ which arise at the two-loop level. Furthermore, we separate the 
non-equilibrium condensate fields: 

a(x) = v + a s {x) 

r(x) = 7r*(x) , z = l,...,iV-l . (5) 

• v is the equilibrium value of the chiral condensate chosen along the sigma direction; 

• a s and 7r* are slow fluctuations about equilibrium. 

The thermal equilibrium ensemble of the hard fluctuations is affected by the presence 
of the condensate. The full, non-equilibrium ensemble averages are basically two- and 
three-point functions of thermalized fields evaluated at the same space-time point. These 
can be written as the sum of the equilibrium ensemble average and a fluctuation about 
this: 



<*?> = + <%/> 

K 2 ) = (vr} 2 ) eg + 5(vr} 2 ) 
((T/tt}) =5(a f n)) 
(ny f ) = 5(n}7f}) 

<*?> = + 5(° 3 f) 

(a fir) 2 ) = {a f u) 2 ) eq + 5{a f Tx) 2 ) 

(Try/) = 5(7r}7rj 2 ) 

(ajvr}) = 5(ay f ) . (6) 

The deviations in the fluctuations (of the general form Sfyfip™)) are the responses of 
the fast modes to the presence of slow a s and tt s background fields. These responses 
are proportional to the slow fields raised to some positive power, and so they vanish, 
as they should, in the absence of the background. Notice the absence of (afH l f) eq and 
(7Tj7Tj) e q. Reason for this is that the correlation functions in equilibrium, in the absence 
of a background, are diagonal. On account of the cubic couplings, however, one finds 
nonzero {cj) eq and ((JfTT^ ) eq . 

In equilibrium, at fixed temperature, a s = and 7? s = and the equilibrium conden- 
sate satisfies the following equation: 

JV-l N-l 

Xv 3 + 3Xv(a 2 f ) eq + XvJ2 (4 2 )«« + X (°f)«i + A E (*f*f 2 U ~ Xv o v ~ H = ■ ( 7 ) 

i=X i=l 
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Since all components of the pion field are equivalent, we can write the sum as X^Ll 1 ( 7r / 2 )eg = 
(N — l)(7Tj) e9 , where now 717 is one of the components. This is a convention that we adopt 
throughout the rest of the paper. 

The resulting final equations describing the evolution of the condensate fields are a set 
of coupled non-linear integro-differential equations. Let us assume only a small deviation 
from equilibrium, |er s |, |7r s | <C v, and thus neglect terms that are higher order in a s and 
-k s . The linearized field equations read 

d 2 a s + M 2 a s + Xv [35(a 2 f ) + (N - 1)8 (rf)] + A \5(a}) + (N — l)5(a f 7r 2 f )] = , 



M 2 = X 2v 2 



f/eq 



(N 



f n f/eq 
V 



H 

+ - , 

v 



(8) 



and 



d 2 n s + M 2 ir s + 2XvS(a f 7T f ) + A 5(ir 3 f ) + 5(a 2 f ir f ) = , 



M 2 w = x(2(7r 2 f ) eq -2(a 2 f ) 



eq 



(9) 



We denote the effective masses of the sigma meson and the pion by M a and M n , and their 
corresponding vacuum value by m a and m n , respectively. 
At one-loop order the effective equations become simpler: 



d 2 a s + M 2 a a + Xv 35(a 2 f ) + (N - l)6(n 2 f ) 







(10) 



M 2 = 2Xv 2 + 



H 



(11) 



and 



Ml 



H 



d 2 n s + M 2 n s + 2XvS(a f n f ) = 



1 ~ 2 /t 2i~2 ~2 o \ T / 2 \ / 2 \ 

+ m n = —m 7r + m n , m n = 2X {7T f ) eq - {a f ) eq 



(12) 
(13) 



and the equilibrium condensate satisfies 

Xv 3 + X \3(a}) eq + (N- l)(ir 2 f ) eq ] v - Xv 2 v -H = 



(14) 
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3 Non-Equilibrium Fluctuations 



Linear response theory is a very convenient tool when one is interested in monitoring 
the effect of an external field applied to a system initially in thermal equilibrium. When 
this effect is small a linear approximation is feasible. A response function expresses the 
difference between the expectation value of an operator before and after an external 
perturbation has been turned on. The effects of the soft condensate fields on the thermal 
medium are thus described by expressions of the form 

8{a%x)ity{x)) = % fdt' f d 3 x'([H perturbed (x'),cr](x)7r] l (x)}} eq . (15) 



o 

Accordingly, we evaluate the commutators of different powers of the field operators at 
two separate space-time points in the fully interacting, but unperturbed ensemble. Initial 
conditions are set by the assumption that in heavy-ion collisions the system reaches a 
state of approximate local thermal equilibrium ||, then it cools while expanding out of 
equilibrium and reaches the critical temperature. We define the initial time t = by 
when the critical temperature is reached. This implies cr s (0, x) — and tt s (0, x) = 0. The 
respective powers n and m are identified from the expressions in equations (§) and (|]). Let 
us emphasize that these responses should be evaluated in the unperturbed, equilibrium 
ensemble, which does include all the interactions between the different modes. 

The possible couplings between low and high frequency modes are determined by 
evaluating the potential U, with the fields separated into their slow and fast components 
as above. The resulting Hamiltonian contains positive powers of the slow fields. For small 
departures from equilibrium it is enough to keep the dominant linear terms only. Relaxing 
the assumption of only slightly out of equilibrium requires the inclusion of higher powers 
of the non-thermal fields. This should provide no difficulties, but is beyond the aim of 
the present paper. The relevant couplings are 

H as a f = A(3w 2 - fl)a s a f + 3\vcr s crj + \a s aj 
H as 7T f = \v<j s ir 2 f 

H-KsVfTTf = A7T s 0j7T/ + 2\viT s a fTYf 

Ha s a f TT f = ^sUfrf (16) 

As before, ir s j refer to one of the N — 1 components of the pion field. 



The response functions obtained by inserting (|16j ) in fll5|) are 



o 



5(aj(x)} = 3i\v / dt' / d 6 x'a s {x'){ aj(x'),aj(x) ) 



ieq. ■> 







5(iTj(x)} = iXv / dt' d x a s (x')( 7Tj(x') } iTf(x) ) eq , , 



o 



S(a f (x)n f (x)) = 2i\v / dt' d 3 x'7r s (x')([a f {x')iT f (x'),cr f (x)TT f (x)}) eq . . (17) 
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The cubic functions are determined analogously, therefore we skip presenting their evalua- 
tion. It is clear that expressions (|17]) vanish in the absence of the non-thermal background, 
as they should. 

The expectation values of the commutators are 



°7 




>°7 


»_ 


)eq. 


II 


ft 




»*/ 


»_ 


)eq. 





D>(x,x f ) 2 , ■ 



([a f (x')n f (x'),a f (x)-K f (x)}) eq . = D<{x, x')D<{x, x') - D> (x,x')D>{x,x') , 

where x = (t,x) and x' = (t',x') are four-vectors in coordinate space. Keeping in mind 
that t' < t, where t is the time elapsed after switching on the perturbation, and t' is 
the time-variable that has its values in the [0, t] interval, the following notation has been 
introduced: 



D>(x,x') = (a f (x)a f (x')) eq . , D>(x, x') = (ir f (x)ir f (x!)) 



eq. 



and 



D<(x,x') = (a f (x')a f (x)) eq . , D<(x, x') = (iT f (x')iT f (x)) eq . 
The functions D > and D K define the spectral function 

p(k) = D > {k)-D < {k) , 
which determines the real time propagator 

d 4 k 



(19) 



D(x,x') = J^e-^-^D(k) 



(2^) 4 
d 4 k 



e -ik(x- X ') (e{t-if) + f(k j)p(k) . 



(20) 



Here k = (k°, k) is the four-momentum, f(k°) = (e fc0/3 — 1) 1 is the Bose-Einstein distri- 
bution and f3 = T~ 4 . The response functions are then 



5(a%x)) 



i6Xv / d x'a s (x') 



d 4 p r d 4 q 



-i(p+q)(x-x') 



(2tt) 4 7 (2tt) 4 

x^(pK(g)(l + /(p°) + /(g )) 



5(tt}(x)) 



i2Xv / d x'a s (x') 



d 4 p r d 4 q 



-i(p+q)(x — x') 



(2tt) 4 7 (2tt) 4 

xp.(p)p.(g)(l + /(p°) + /(g )) 



5{af(x)irf(x)) = —i2Xv / d x'7r s (x') 



d A p f d 4 q 



-i(p+q)(x-x') 



(2tt) 4 7 (2tt) 4 

xpAp)Pn(q)a + f(p°) + f(q )) ■ (2i) 



9 



4 The Order Parameter and Meson Masses 



To analyze the phase transition we look at the temperature dependence of the equilibrium 
condensate which is the order parameter of our model. The behavior of v is determined 
by solving equation (|7|). It is educational to look at the theory in the exact chiral limit, 
H = 0, first. In this case equation (0) for the order parameter has two solutions: 

v = , T > T c , 

and 



v 2 



fl ~ ZtfU -(N- l)(nj) eq - 'flhl - (jV - l)^hl , T < T c . (22) 

v v 



The first solution shows, as we expect, that the chiral condensate does not exist above 
a critical temperature. The second solution represents the low temperature, symmetry 
broken phase. The two solutions are separated by a second order phase transition. 

The order parameter depends on the equilibrium field fluctuations, which to first 
approximation in either perturbative expansion in A |yj or 1/iV-expansion [27[], are 



<^>- = /(&^ 1+2 « B »' (23) 

where (p = 07,717 and / is the Bose-Einstein distribution function of the mesons with 
energy E = \fp 2 + m 2 . Cubic fluctuation terms that enter in the expression for the 
equilibrium condensate, (a^) eq /v and {afn 2 ) eq /v are non-zero on account of the possible 
cubic couplings. However, already the leading order term of the expansion is higher order 
in the coupling than (^), so we can drop these to first order. 

Let us discuss some important observations regarding the fluctuations (^). First, 
keeping only the leading order term means that the masses are the bare zero temperature 
masses. However, as indicated by equations (|TTJ) and flTB|), the masses themselves are 
given in terms of the equilibrium condensate and as such, they are temperature dependent. 
Therefore, the temperature dependence of the order parameter should be determined by 
performing a self-consistent evaluation. Here we resum tadpole diagrams. 

Second, the first term is the vacuum contribution, while the second term is due to 
finite temperature effects. The zero temperature part is divergent in the ultraviolet limit. 
This divergence can and should be removed using vacuum renormalization techniques. 
The finite temperature does not introduce any extra divergence since it is regularized 
by the distribution function. One can say that the very short distance behavior of the 



theory is not affected by finite temperature |3T[]. Therefore, T = renormalization is 
enough to obtain finite results. The usual approach then is to neglect the zero temper- 
ature contributions when focusing on the physics at finite temperatures. At this point, 
it is worth mentioning that using a self-consistent approximation makes the usual renor- 



malization procedure difficult, as discussed in ||32|| . The argument, according to which the 
renormalized divergent term can be neglected, is really correct only when the mass is the 
bare mass. A self-consistent calculation involves the temperature-dependent mass, lead- 
ing to the temperature-dependence of the divergent term. Renormalization thus results 
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in temperature-dependent renormalization constants, and these should not be ignored. 
However, such a treatment is beyond the scope of this paper, and in what follows, we are 
going to ignore the divergent term, while still being aware of this approximation. 

Third, the momentum integration has a lower limit, A c , due to the restriction of 
the thermal population to hard momenta, \p\ > A c . When A c = 0, the condensate 
contains only zero momentum modes, meaning that the classical field configurations are 
homogeneous. In reality non-zero but small momenta can be part of the condensate. 
Then we talk about inhomogeneous condensate. 



4.1 Goldstone Modes 

In the theory with spontaneously broken chiral symmetry the tree level mass of the pion 



is zero in the broken phase, m % = 0. Goldstone's Theorem |23| requires that this remains 
zero at every order in perturbation theory. The first glimpse at equation (H) (or (|13"D ) for 
the pion condensate shows a mass-term. At one-loop order 



ml 



2A 



<*?> 



f/eg 



<*?> 



f/eg 



7^0 



(24) 



which is zero only at i) zero temperature, where the thermal fluctuations themselves 
vanish, or ii) when the masses of the pion and the sigma are equal, which is expected at 
the critical temperature, or iii) in the high temperature limit, when the meson masses 
can be neglected. In the following, we present a simple and clear way to prove that this 
violation of Goldstone's Theorem is only apparent. 

The pion mass ( p4|) includes one-loop tadpole contributions. There is another one- 
loop diagram that contributes to order A, 'built' out of two 3- vertices, known as sunset 
diagram. This diagram is incorporated in the equation of motion through the response 
function 5(iVf(Tf) given by fl2"ip. The function is of the order of Xv ~ A 1//2 , and there is 
an overall factor of 2Xv ~ A 1 / 2 in front of it in the effective equation. Therefore, to get 
the contribution of order A, the expectation value in (|2l|) can be evaluated at the lowest 
order. This means that we can replace the interacting ensemble average (...) e <j. by the free 
ensemble average (...)o- This in turn is equivalent to inserting free spectral functions, 



Pfree{p) = 27ie(p°)5 (p 2 - 771 2 ) , 

into the expressions of fl21~|) . After evaluating the frequency integrals we find 



5(a f n f ) = ,2AM 2 / d'x'^x') J ** j 



d 3 p f d 3 q i 



(p+q)(x-x ') 



(2tt) 3 7 (2tt) 3 



(25) 



(26) 



with 



F(p,q,t') 



1 



[1 + U + U) {e l(E ° 

+(/* - U) ( 



+E % )(t-t') 



i{E„-E v ){t-1/) _ -i(E a -E w )(t-t) 
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Since, the deviation from equilibrium is assumed to be small, we Taylor expand k s {x') 
about its equilibrium value. The first term of the expansion, linear in k s (x), gives the 
contribution to the mass. The total pion mass is then 



M 2 = ml + cii 

7T 7T L 



where 



«i 



i2\M 2 a 
2\ 



d 3 p 
(2tt) 3 
d 3 p 



dt'F(p, t') 



1 



d 3 p 



(2tt) 3 4E a E v 
1 



2E a {l + 2f v ) - 2E v (l + 2j a ) 



(2k) 



2E n 



2E n 



m - ED 



and E a = ^Jp 2 + Ml and ^ 
true pion mass yields 



p |. Summing up ( p8f ) and in expression 



M T = 



(27) 



(28) 
, the 

(29) 



thus proving that with proper inclusion of diagrams the pions stay Goldstone bosons at 
one-loop level at all temperatures. Another simple and straightforward proof in frequency- 
momentum space is presented in p2|. 



4.2 Numerical Results 

In the exact chiral limit the self-consistent solution of the gap equation 



v 



2 --fl-^a}) eq -(N-l)^)) eq , (30) 



with the sigma and pion field fluctuations 

1 r°° 

7 ' ) = 2^7a c aiJ E~e E °/ T -l 



( a f) = TT"? / d P 7^ W~Tt 7 > ( 31 ) 



and 

1 r°° v 2 1 T 2 
^ > = 2^1*fe^7^T^12- fi»M.^0,A.-0, (32) 



evaluated with E a = y p 2 + Ml, where 

Ml = 2Xv 2 , (33) 

is presented in figure [l] for different values of the momentum separation scale A c . In 
the numerical analysis N = 4 and the coupling constant was chosen to be A = 18, 
corresponding to a vacuum sigma mass of about m 2 = 2Xfl = (558 MeV) 2 . The phase 
transition temperature, defined by the vanishing of the condensate, is about T c ~ 130 MeV 
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Figure 1: Temperature dependence of the sigma mass for different momentum separation 
scales in the chiral limit. 
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Figure 2: Temperature dependence of the equilibrium condensate for different coupling 
constants in the chiral limit. 
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for homogeneous classical field configurations, A c = 0. With the increase of the scale A c 
the phase-space available for hard modes is decreased, requiring higher temperatures for 
the fluctuations to completely dissolve the condensate. For A c = 120 MeV the temperature 
above which fluctuations are much too large to allow the formation of any condensate is 

Figure |I| shows 



about T c ~ 



165 MeV, which is in good agreement with lattice data [33 



that the meson mass is positive-definite at all temperatures, eliminating the tachyon 
problem present in the mean field approximations [[31] . 




Figure 3: Temperature dependence of the resummed meson masses, and of the equilibrium 
condensate in the 0(4) model with explicitly broken symmetry. 

The theory with exact chiral symmetry is known to have a second order phase tran- 
sition based on universality arguments pi. This has been confirmed within the mean 



field approximation ||27|| . Figure [I] however, shows a discontinuous behavior of the order 
parameter at T c . Such a jump is characteristic to first-order phase transitions. Incorpo- 
rating the effect of thermal fluctuations in a self-consistent way renders the transition first 
order. Such a behavior has been discussed many years ago also by Baym and Grinstein 
26j . In figure |2] we present how the discontinuity decreases with decreasing coupling con- 



stant. For a weak enough coupling the continuous second-order transition is recovered. 
This can be understood in the following way: in order to assure the finiteness of the 
0(iV)-theory in the large- N limit the coupling constant should be written as X/N |26j in 
the self-consistent gap equation fl3"5|)-(|3"(]|): 



Ml = m 2 -2- 
N 



3(a 2 f ) eq + (N-l)(rr 2 f ) eq . (34) 



For N — > oo the contribution from the sigma field fluctuation disappears. The condensate 
is then 



V^fl-^WU^fl-^ (35) 
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clearly showing a continuous decrease of the order parameter with increasing temperature. 
In the 0(4) model decreasing A by hand is equivalent to going to the large N limit in 
the 0(N) model. For our model parameters we have determined a large N CT iu C ai — 1800 
at which the transition is second order. For this A = 18 was held fixed. This result is 
equivalent to having X cr iucai — 0-01 and N = 4. 

Numerically determined self-consistent solutions for the condensate ( II ) and the me- 
son masses (|IID and (|T3|) in the more realistic theory with explicitly broken symmetry 
are displayed in figure |3|, and show a qualitatively different behavior. There is no phase 
transition in the textbook sense. The equilibrium condensate monotonically decreases 
with increasing temperature. A crossover region can be defined where the sigma and pion 
masses start to approach degeneracy. The minimum of the sigma mass is at T ~ 235 MeV. 
Different values for A c do not introduce significant effect in evaluating the meson masses. 



5 Dissipation of the Chiral Condensate 

Dissipation of the condensate occurs because energy can be transferred between the con- 
densate and the heat bath through the interactions of soft and hard degrees of freedom. 
Formally, in our model, the damping of different modes is determined from the response 
functions. Most of previous studies on this topic have been done at one-loop level. The 
two-loop level scattering processes have been ignored, based on the argument that these 
are higher order in the coupling constant than are the decay and absorbtion processes. 
However, as we show in the following, two-loop processes can dominate due to the available 
large phase space. 

The effective equations of motion for long-wavelength meson fields (jl(|) and (|T2"D have 
the general form 

d 2 <p s (x) + M 2 <p s (x) + F(x) = , (36) 

where S = a s or 7r s and 

F(x) = J d*x'(f) s (x')Tl(x -x') . (37) 
In frequency-momentum space this reads as 

— k 2 + M 2 + n(Jfe) = . (38) 

The function Tl(k) is given by 

n(jfe) = -ig 2 I d 4 x' [ [ _^ e <(*-p-*)(*-*') 

K J y J J (2vr) 4 i (2vr) 4 

xpi(p)p 2 (g)(l + /(p°) + /(g°)) , (39) 

and is identified as the self-energy. Here g is the corresponding coupling, and the indices 
1 and 2 refer to either of the hard modes, aj, and iif, respectively, and k = (k°,k) is 
the four-momentum of the soft sigma meson or pion. The frequency has a real and an 
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imaginary part, k° = u ■ 
in the dispersion relation 



iT, provided k real. The real part of the self-energy participates 



k 2 + M 2 + Re IL(u, k) 



(40) 



and with the usual assumption of weak damping, r C w, the imaginary part of the 
self-energy completely determines the damping of excitations: 



r 



Im H(uj, k) 
2oj 



(41) 



T is the rate at which out-of-equilibrium meson modes with energy u and momentum 
k approach equilibrium. This rate is determined by physical processes which can be 
identified from the imaginary part of the self-energy ||34|| . One should be aware that (|4l| ) 
describes the rate of decay of the amplitude of the wave, exp(— Tt). The loss rate for the 
number density is exp(— 2Tt). 



5.1 Dissipation at One-Loop Order 

There are several diagrams contributing to the self-energy at one-loop order. Tadpoles 
are real and they only modify the mass, as we discussed before. Dissipative effects come 
from non-local diagrams. At order A these are determined from the response functions 
through ( PU|) with the insertion of the free spectral functions, 

d 3 p 1 



U(k) =g 2 J 



(2n) 3 AE 1 E 2 

+(/ 2 



1 + fx + h 
fx 



1 1 

uj — Ei — E 2 + ie uj + Ei + E 2 + ie 
1 1 



uj — Ei + E 2 + ie uj + Ei — E 2 + ie 



with Ei = ym 2 + (p + k) 2 and Ei = \Jm 2 + p 2 . The indices 1 and 2 refer again to either 
of the fast sigma and pion. It is important to observe that the above expression coincides 
with the self-energy calculated directly from the non-local one-loop diagram [33]. For 
positive energies of the mesons u > the contributions to the imaginary part of the 
self-energy are 

d 3 p 



Im Tl(uj, k) = —Tig' 



1 



{2ufAEiE 2 



[(1 + fx + / 2 )% -Ex- E 2 ) 

+(f2-fx)5(uj-Ex + E 2 )] 



(42) 



Because the heat bath singles out a preferred reference frame we keep u and k independent. 
The dynamics of the decay is determined by the on-shell processes that are allowed. 



Sigma Meson Decay. Contribution to the imaginary part of the sigma meson self- 
energy comes from the diagram presented in figure Kinematically, the decay of a soft 
sigma meson into hard thermal pions, a s — > ir fir f, and the inverse, recombination process 
is allowed, provided that uj 2 — k 2 > Am 2 . 



Im H(uj, k) 



3_ 

IQir 




(43) 



16 



(J, V 



(7c 



Figure 4: One-loop self-energy contribution to the soft <r s with coupling to the hard ttj through 
Xva s ir 2 f . 



where 




1 - 



4m 2 



k 2 



The decay and formation processes are obtained by setting u 2 — k 2 



mz 



Im U(m n 



167T 



\ 



( 



1 - 



4m 2 7 1 - exp v 2;r 2 , 



+ 2- log 



/ 1 \ \ 



\ 



exp(-^ + ^Wl - 



(44) 



The rate at which soft, k <T, sigmas of energy u; disappear from the condensate due to 
their decay into thermal pions is 



a; 



(N — 1) ml — ml 



16n 



-X- 



A 



l^cothf^ 
m 2 V4T 



(45) 



The temperature dependence of the sigma damping rate (^) in the rest frame of the 
is shown in figure [5| The calculations were done with the resummed 



sigma, u = m 

meson masses, M a and M n . Figure [5] shows that even at zero temperature there is a 
finite damping, so the sigma meson can decay into two pions in the vacuum. At T = 0, 
for our model parameters the damping rate is about T a7T7T = 510 MeV, which is of the 
order of the mass of the sigma. The width of the sigma resonance is very broad, in other 
words the sigma meson is overdamped. Figure |5| shows that the damping is increasing 
with T, and is followed by a quite drastic decrease starting from about T = 150 MeV. 
This temperature corresponds to the temperature where the sigma mass begins to drop 
significantly (see figure |3|). It is then natural to expect a decrease of the sigma decay rate 
into pions that have masses approaching that of the sigma. Above T ~ 200 MeV the 
threshold condition for the decay of an on-shell sigma meson, M a > 2M n , is not fulfilled 
anymore. Therefore, there will be no contribution from decay to the sigma width in the 
kinematically suppressed region. 
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Figure 5: Temperature dependence of the sigma damping rate at one- loop order in the 
sigma rest frame. Calculations were done with resummed meson masses. 



Pion Damping. At one-loop order there is only one diagram contributing to the pion 
self-energy. This diagram is shown in figure |6|. At this order, dissipation of the pion 
condensate can occur provided the energy and momentum of the soft pion satisfies the 
kinematic condition uo 2 — k 2 < (rricr — m^) 2 . Then the transformation of a pion into a sigma 
when propagating through a thermal medium can happen. Basically, a soft pion from the 




Figure 6: One-loop self-energy contribution to ir s with coupling to a s and tt s through 2\vir s afKf. 



condensate annihilates with a hard thermal pion producing a hard thermal sigma meson, 
7T s iTf — > 07. The inverse process is the decay of a hard thermal sigma meson into a soft 
and a hard pion, <jf — > 7i s 7if. The net rate of dissipation is 



8ir ku 



UJ I + UJ 



log 



e t 



e t — 1 



+ log 



e t 



e t 



(46) 



where 



p A ± + m% 
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p± 



k m 2 — 2m 2 



mi 



u> m a 
2ml\ 



Ami 



mi 



In the rest-frame of the massive pion the expression for the damping is simplified to 



A ml(m 2 — m 2 ) 



16tt 



mi 



A 



4ml 



T 



mi 



""a 

1)(1 - e"2^r) 



(47) 



Figure [F] shows the temperature dependence of the damping rate of massive pions at one- 
loop order fl4"?D calculated in the pion's rest frame using the resummed meson masses, 
M a and M n . Note that at zero temperature the dissipation is zero. This makes sense 
because the transformation of pions into sigmas is due to their annihilation with a hard 
thermal pion in the medium, and so this is exclusively a finite temperature process. At low 




Figure 7: Temperature dependence of the pion damping rate at one-loop order in the pion 
rest frame. Calculations were done with resummed meson masses. 



temperatures the phase space available for this pion transformation process is suppressed 
by the large sigma mass. As the temperature increases and the sigma mass is dropping 
the width of the pions is increasing. Figure [7| shows that the damping can get quite 
strong. At T ~ 170 MeV, for example, when the pion mass is about m n = 158 MeV the 
damping is T nna — 87.0 MeV. Thus, at this temperature the width of the pion is about 
55% of its energy. This result makes us question whether the pion is a quasiparticle in this 
temperature region and needs further investigations. At temperatures around 200 MeV 
the kinematic condition for an on-shell pion, M a > 2M n , is not satisfied, prohibiting the 
transformation of a pion into a sigma while passing through a hot medium. 

The damping in terms of the pion momenta is presented in figure ||] for different tem- 
peratures. Results using the dispersion relation u 2 = k 2 + M 2 are displayed. Noticeable 
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Figure 8: Momentum dependence of the pion energy and width at different temperatures 
calculated with resummed meson masses and with dispersion relation uj 2 = k 2 + M 2 . All 
quantities are normalized to the vacuum pion mass. 



damping occurs above T = 100 MeV and increases with T. At about T = 160 MeV all 
modes are equally damped. In other words, the width of the pion is independent of its 
momentum. This width is increasing with T and it can be as great as 30% of the energy. 
At even higher temperatures the damping of the zero momentum modes is the strongest. 

5.2 Dissipation at Two-Loop Order 

Instead of the tedious evaluation of two-loop linear response functions, we present a 
direct determination of damping rates from the physical processes responsible. Two- 
loop contributions correspond to two-particle scatterings with amplitudes evaluated at 
tree level. Similar to the one-loop calculations we evaluate the imaginary parts of the 



two-loop self-energies and inserting these into (|4T|) results in the damping rate due to 
scattering. 

The general form of the self-energy of a particle of mass m a , propagating with four- 
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momentum k = (u, k) through a medium in thermal equilibrium, is given by 

^'^/pfb^w- (48 > 

Here Ai is the transition amplitude for the scattering process ab — > ab. The ther- 
modynamical weight f(E) is the Bose distribution of thermal mesons of mass m b and 
four- momentum p = (E,p). In terms of the forward scattering amplitude Ai(s) = 
— 87r v /i/ cm (s), where s = (p + k) 2 , and the imaginary part follows 

lmU ab (k) = -J^- 3 f(E)V~s q -^a total (s). (49) 

To obtain this equality we applied the standard form of the optical theorem that relates 
the imaginary part of the forward scattering amplitude and the total cross-section [36]: 

Im f cm {s) = -pC7totaf(s) • (50) 

Here we consider scatterings involving massive particles only. Then 

Im U ab (u) = - JL J' cos 9 J™ dp^f(E)a ab (E)^(s - m g + m 2)2 _ 4sm 2 ? 

with s = m 2 a + m 2 + 2Eu — 2pkcos8, where 9 is the angle between k and p. The 



dispersion relation of the hard thermal modes is E = Jp 2 + m 2 and of the mesons inside 



the condensate uo = yk 2 + m 2 . It is convenient to work in the rest frame of meson a, 
where 

irn poo 

Im U ab (uj = m a , k = 0) = --| / dE(E 2 - m 2 b )f(E)a ab (E) . (51) 
The cross section for a scattering ab — > a& is given by 

\ / cm 

where 

/rf(TaA 1 \M\ 2 . (53) 



The symmetry factor 1/ S\ is due to the number S of identical particles in the final state. 
It is clear that knowing the amplitude M. of a process readily results in the dissipation 
rate due to that process. 
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Sigma Scattering. There are two mechanisms that contribute to the removal or addi- 
tion of a sigma meson to the condensate: elastic scattering of a hard thermal sigma or 
of a hard thermal pion off a sigma from the condensate. To first order in the coupling 
A there are four diagrams contributing to the process in which a thermal sigma meson 
knocks out a low momentum sigma from the condensate. The transition amplitude, the 
sum of contributions from different diagrams is: 



M = -6A 



1 + 3(m 2 — m 



1 



t 



+ 



u 



m 



(54) 



This reflects the symmetry in the s, t, and u channels, and s + 1 + u = Am 2 a . The total 
cross-section is 



On 



is) 



9A 2 

8ns 



' s + 2m 2 -3m 2 



2\2 



+ 



mi 



18(m 2 -m 2 ) 
m 2 (s — 3m 2 ) 



12(m 2 - m 2 )(s 2 - 3sm 2 - m% + 3m 2 ml) 



In 



3m 2 



(s — Am 2 )(s — 2m 2 )(s — m 2 ) 
In the low energy limit expand this about s = <Am 2 a . To leading order 



(55) 
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32~k' 



A 2 



(4m 2 — 5m 2 ) 2 



mi 



and in the rest-frame of the sigma this gives rise to 



Im IL 



9 



2\2 



A 2 T 2 e -m CT 



327T 3 



/T (m a + T)(4m a -5m 7T ) 



mi 



(56) 



(57) 



In the high energy limit only the four-point vertex contributes to the amplitude, resulting 
in 

9A 2 



8ns 



This gives rise to 



Im IL 



2rji2 



64tt 



AT 



(58) 



(59) 



With these two limits we construct an interpolating formula that describes the whole 
energy range. The contribution to the rate of decay of the amplitude is then 



r ~ 



9A 2 



T 2 (m CT + T)(4m 2 -5m 2 



2\2 



64tt 6m a (m a + T)(4m2 - 5m2) 2 + ix 2 m%{e m °l T - 1) 



(60) 



A hard pion in the heat bath can be energetic enough to knock out a sigma meson 
from the condensate. The possible tree-level processes may happen according to the four 
different diagrams. The transition amplitude obtained from these is 



M = -2A 



1 + (m 2 - m 2 n ) 



+ 



+ 



mi 



mi 



u 



mi 



(61) 
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where s + t + u = 2(m 2 + m 2 ) . The total scattering cross-section is 



On 



4-7TS 



s — 2m 2 + 771 ^- 



9s(m 2 — ml) 2 



mt 



s[m„ 



ml(s 2 — sm 2 a — 2sm 2 + (m 2 — m 



2 \2>, 



2\2 



(2m 2 + m 2 — s)((m 2 — ml) 



smt 



6s(m 2 — m 2 )(sm 2 + 2sml — s 2 + — — 2m 2 m 2 ) 
(s — m 2 )(m 2 + ml — s)(s 2 — 2s(m 2 + m 2 ) + (m 2 — ml) 2 ) 

X l n ^1 

s 2 — sm 2 — 2sm 2 + (m 2 — m 2 ) 2 

2s(m 2 — m 2 )(3sm 2 — s 2 + + m\ — 4m 2 m 2 ) 
(s — m 2 )(m 2 + ml — s)(s 2 — 2s(m 2 + m 2 ) + (m 2 — m 2 ) 2 ) 
s(2m 2 + m 2 - s) 



x In 



fm 2 — m 2 ) 2 — sm 2 



(62) 



The low energy limit is obtained by expanding this about s = [m a + m^) 2 and is 

= 9 A2 m^(3m 2 -4m 2 ) 2 ^ 
47r sm 2 (m 2 — 4m 2 .) 2 



(63) 



At high temperatures, where only the four-vertex diagram contributes, the cross-section 
is reduced to 



<7n 



4tTS 



The contribution to the imaginary part of the self-energy at low energies is 

9 



Im IL 



._ X 2 T 2 e -m n /T F ^ m ^ m ^ ; 



where we defined 



F U7T (m a ,m n ) 



2\2 



ml(m n + T)(3m 2 — 4m 2 
m z a {m a + m 7r ) 2 (m 2 — 4ml) 2 



and at high energies 



\ 2 T 2 
96vr 



im n r 



(64) 



(65) 



(66) 



(67) 



The two limits can be combined into one approximate expression which then determines 
the rate of dissipation. Factoring in all the N — 1 pions 



^ 9(N — 1) 2 T 2 

-L /T7T „ 'V 



S7T 



(6? 



m a 216F a7T + 7r 2 (e m -/ T - 1) 

Distinction should be made between scatterings with massless and massive pions. The 
low energy expression (53) vanishes for zero pion mass. The first nonzero term in the 
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series expansion of the cross section is the fourth order term. We discuss the massless 
case in our forthcoming paper p2 . 

The temperature dependence of the total scattering rate of the sigma meson, T aa + 
T a7T , evaluated with the self-consistently determined meson masses is shown in figure []. 
Scattering is more accentuated at higher temperatures and its contribution to the sigma 
damping rates is well below the energy. Dissipation of sigmas from the condensate due 
to their scattering is much smaller than due to their decay, meaning that the scalar order 
parameter relaxes to its equilibrium value via the production of lighter pion fields. It also 
means that sigma mesons are so unstable that they are more likely to decay before they 
could ever scatter with other particles from the medium. 
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Figure 9: Scattering contribution to the width of the sigma meson as function of temper- 
ature. Calculations were done with resummed meson masses. 



Pion Scattering. Dissipation of DCCs can arise from scattering the of soft pions with 
hard pions or hard sigma mesons. The damping of massive pions and that of the Goldstone 
pions we expect to be different, requiring a somewhat different analysis. In the following 
we present the discussion on massive pions. There are four possible tree-level diagrams 
representing the reaction in which a hard thermal sigma knocks out a low momentum 
pion from the condensate. The total cross section is the same as for sigma-pion scattering 
and is given by expression (|62|) . The imaginary part of the self-energy in the low energy 
limit in the rest frame of the pion is 

Im n™ = - A A 2 TV m ^ T F™ , (69) 



where 



= ml{m a + T){?,ml-Aml) 2 
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and for high energies is 



A 2 T 2 
96vr 



im n T 



(71) 



The total scattering rate due to this process can be parametrized by an interpolating 
formula between the two known limits: 



9A 2 T 2 



r ~ 

- 1 - 7T<X 



87rm ff 216F 7Ta + n 2 (e m °/ T - 1) 



(72) 



im scattering has been extensively studied during the last couple of decades in a 
variety of different models and approaches. An incomplete but significant list of references 
is |37]. Here we study the elastic scattering of a hard pion off a soft pion. For one pion 



species, N = 2, there are four diagrams contributing to this process at tree-level: one 
4-point vertex diagram and three 3-point vertex contributions involving a sigma exchange 
in the s, t and u channels. The transition amplitude is 
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(73) 



where s + 1 + u = Am 2 . Accounting for all pion species opens up additional channels. For 
N = 4 the transition amplitude averaged over initial and summed over final isospins is 



\M\ 2 = \ £ (2/+l)|A^ 2 

6 1=0,1,2 



(74) 



where M 1 is the matrix element associated with the total isospin I of the two-pion system: 
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The resulting cross-section is 
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(75) 
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The low-energy limit of the cross-section, quite acceptable for m a 3> T, is given by the 
expansion of this about s = 4m 2 , 
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At high temperatures, in the T — > T c limit the major contribution to the amplitude is 
from the four-point vertices, resulting in 

ISA 2 

°w = s — • (< 8 ) 
8ns 

In the rest frame of the pion the above two limits give 

32tt 3 (m 2 -4m 2 ) 2 1 } 

and 

5A 2 T 2 

Im n„ = -— — , (80) 

respectively. The rate of dissipation due to massive pion-pion scattering is then given by 
the interpolating expression 

r ., 23A 2 T 2 ml(m n + T) 

64tt f m3(m^ + T) + vr 2 (m 2 - 4m 2 ) 2 (e^/ T - 1) ' 1 ' 




Figure 10: Scattering contribution to the width of the pion as function of temperature. 
Calculations were done with resummed meson masses. 

The scattering contribution to the damping of massive pions is presented in figure 



roj . Due to the heavy sigma exchange there is a strong suppression at low temperatures. 
When reaching T ~ 130 MeV the scattering rate becomes significant and comparable 
in magnitude to the damping at one-loop order. In the critical region the contribution 
to the pion width from pion-pion scattering grows rapidly reaching a maximum about 
T = 200 MeV. It is interesting to note from figures ^| and [Tt] that in the high temperature 
regions, where the mesons become almost degenerate their widths approach the same 
value, as expected. 
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6 Relaxation Time 



The ability to detect DCCs in relativistic collisions of heavy ions depends on the lifetime 
of the condensate. DCC formation can happen out of thermal equilibrium only In order 
to talk about non-equilibrium physics the rate of expansion,t ea; p, has to be much smaller 
than the relaxation time of long-wavelength modes, t, 

texp "C t . (82) 

Otherwise, for a slow expansion, the soft modes have enough time to equilibrate. Equili- 
bration of the out-of-equilibrium chiral condensate is the result of the presence of a heat 
bath. Above, we analyzed the physical processes responsible, and calculated the damp- 
ing of different meson modes. The total width is the sum of one- and two-loop order 
dissipative contributions, 

l-\r ~^TTn<7 ~t~ 1^770" ~t~ f^TT 5 (^^) 

which exhibits a sharp peak in the critical region, due to the peak in the damping rate 
from pion-pion scattering. In the non-equilibrium, but close to equilibrium physics that 
we consider, the damping directly controls the rate at which equilibrium is approached 
through the relation [|34| 

* = ^ ■ (84) 
Obviously, the larger the damping due to the interaction of the condensate with the heat 



bath, the shorter the relaxation time is. Figure [LI] shows the change in the relaxation time 



with temperature. At low temperatures t is large due to the suppression of the thermal 
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Figure 11: Relaxation time of a homogeneous disoriented chiral condensate versus tem- 
perature. 
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occupation numbers. As the temperature increases more thermal modes get excited and 
as a consequence the relaxation time is quickly decreasing. In the phase transition region 
the decay time is the shortest. At the peak of the damping T = 200 MeV we found 
t = 0.17 fm/c, and at T = 235 MeV, where the mass of the sigma is the smallest, 
t = 1.36 fm/c. When assuming no multiple interactions with the heat bath then t is the 
lifetime of the DCC. The time we obtained are shorter than previous estimates [^|, IB], 15], 
and are short enough to make a possible DCC signal questionable. One can expect 
that any multiple scatterings or decays would only increase the damping, decreasing the 
relaxation time. 



7 Conclusions 

This work was motivated by our interest to determine the possibility of survival of DCCs 
in the background of a multitude of thermal particles, mostly pions, that are formed 
after two heavy ions are collided at ultrarelativistic energies. Also, having a consistent 
description of quantum fields near thermal equilibrium allows for a better understanding 
of their dynamics in a region where non-perturbative analysis is required. 

We have developed a consistent semi-classical study of the out-of-equilibrium chiral 
condensate fields in the framework of the linear sigma model. Clear distinction between 
the soft non-thermal chiral fields and hard thermal modes has been made, accounting also 
for interactions between these. Motivated by the high occupancy of the low momentum 
modes we allowed for their classical treatment. The effect of the other degrees of freedom 
has been taken into account by introducing a heat bath of mesons. These thermalised 
high momentum modes have been accounted for in a perturbative manner, improved by 
the resummation of certain diagrams. We derived classical equations of motion for the 
long-wavelength condensate fields coupled to the thermal bath. After integrating out the 
hard modes effective field equations resulted, which completely determine the evolution 
of the chiral condensate in space and time. 

The presence of the slowly varying condensate fields cause deviations in the equilib- 
rium fluctuations of the thermal fields. We identified these as linear response functions, 
since we are dealing with not too far from equilibrium scenarios. We have discussed 
in details the richness of information contained within these response functions: They 
renormalize the equations of motion, modifying the particle properties, and give rise to 
dissipation. 

The temperature dependence of the meson masses and that of the equilibrium con- 
densate have been determined numerically in a self-consistent manner, both in the chiral 
limit and for explicitly broken chiral symmetry. In the chiral limit a first-order phase 
transition was found, which is an artifact of the model. We show that for a small coupling 
constant, or large N limit, the expected second-order transition is recovered. We have 
considered in some detail the Goldstone boson nature of the pion, proving that when 
properly accounting for the tadpole and sunset diagrams the pion remains massless at 
one-loop level in the symmetry broken phase. Also, the tachyon problem of mean-field 
approximations is eliminated by naturally assuring the positivity of the masses at all tem- 
peratures. In the more realistic case, in which chiral symmetry is explicitly broken by the 
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nonzero quark mass, a crossover region was identified. The minimum of the sigma mass 
is at the transition temperature of about 235 MeV. Above this the masses of the pion and 
the sigma become degenerate and the equilibrium condensate vanishes asymptotically, 
signalling an approximate restoration of chiral symmetry. 

Due to possible interactions between different degrees of freedom, those of the conden- 
sate and those of the heat bath, energy exchange is possible and particles can be knocked 
out or put in the condensate. Direct evaluation of the response functions results in the 
rates for the different processes. We have identified these physical processes that are re- 
sponsible for the dissipation of long-wavelength modes of the chiral condensate, and have 
confirmed that at high temperatures not only the damping of the sigmas is significant, 
but also that of the pions. At one-loop level, provided a kinematic condition is satisfied, 
a pion from the condensate can annihilate with a pion from the heat bath forming a 
sigma. The damping due to this process becomes stronger with increasing temperature. 
The width of the pion can be as big as 55% of its energy. This result makes us question 
whether one can talk about the pion as a quasiparticle in the phase transition region. 

We emphasize the importance of two-loop calculations. Our results show that contri- 
bution to the damping rate of the condensate provided by two-particle elastic scattering 
processes can get as important as are decay processes. Moreover, while decays happen 
only when certain kinematic conditions are satisfied elastic scatterings have no such re- 
strictions. Therefore two-loop processes contribute to the width even in regions where 
one-loop contributions are suppressed. Soft pions from the condensate can be knocked 
out through elastic scatterings with a hard thermal pion or sigma. The damping due to 
these scatterings is most accentuated in the phase transition region, and is the greatest 
at about 200 MeV, where the pion width shows a peak. 

The damping directly controls the rate at which equilibrium is achieved by the non- 
equilibrium condensate. We have determined the relaxation time of a homogeneous con- 
densate due to both one- and two-loop order dissipative processes. We have obtained 
relaxation times between 0.17 — 1.36 fm/c in the phase transition region. Assuming no 
multiple interactions these times become the lifetime of the condensate. We have found 
that the lifetime of disoriented chiral condensates is short enough to make a possible DCC 
signal questionable. 

The natural next step of our investigation is to solve the field equations for some 
initial conditions for an expanding system. Such analysis is on the way. 

The methods used in this paper are general, and may be used in other contexts, where 
non-equilibrium physics of quantum fields is of interest. 
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